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ABSTRACT 


Thermal buckling characteristics of hypersonic aircraft sandwich panels of various aspect ratios were 
investigated. The panel is fastened at its four edges to the substructures under four different edge condi¬ 
tions and is subjected to uniform temperature loading. Minimum potential energy theory and finite element 
methods were used to calculate the panel buckling temperatures. The two methods gave fairly close buck¬ 
ling temperatures. However, the finite element method gave slightly lower buckling temperatures than 
those given by the minimum potential energy theory. The reasons for this slight discrepancy in eigensolu- 
tions are discussed in detail. In addition, the effect of eigenshifting on the eigenvalue convergence rate is 
discussed. 
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a x , a y coefficients of thermal expansion, in/in-°F 

a xy coefficients of thermal shear distortion, in/in-°F 

y xz , y yz transverse shear strain in xz- and yz-planes, in/in 

AT temperature rise, °F 

A T a assumed buckling temperature, °F 

A T cr critical buckling temperature, °F 

X, numerical coefficient of N T y in a^ nk[ 

r] numerical factor in buckling equation, which changes with the edge condition 

X i eigenvalue of i-th buckling mode 

v xy , v yz Poisson ratios of face sheets, also used for those of sandwich panel 

v cxy’ v cyz> v cxx Poisson ratios of honeycomb core 

X, numerical coefficient of N x in a^ nkl 

p Ti specific weight of titanium material, lb/in 3 

p Hc specific weight of titanium honeycomb core, lb/in 3 

INTRODUCTION 

Hypersonic aircraft structural panels are subjected not only to aerodynamic loading (mechanical load¬ 
ing), but also to aerodynamic heating (thermal loading). These structural panels are usually called hot 
structural panels because they operate at elevated temperatures. For certain cases, the thermal load could 
be the primary load, and therefore, it could be a key factor in the design of the hot structures. When a 
monolithic hot structure is subjected to uniform temperature field and is allowed to expand freely, no ther¬ 
mal stresses can be generated in the structural panel. When the temperature field is nonuniform, thermal 
stresses can build up in the panel even if it can expand freely. In actual applications, the structural panels 
are attached to relatively cooler substructures (i.e., spars and ribs, both of which act as heat sinks); the pan¬ 
els are, therefore, constrained from free expansion. These constraints will cause thermal stresses to build 
up in the panels. The heating over an individual panel surface is usually relatively uniform; however, the 
panel surface temperatures are seldom uniform over the entire panel surface because the panel edges are 
attached to the heat sinks (i.e., relatively cooler substructures). The temperature rise over the panel surface 
will then have a plateau in a large central region and will taper down near the cooler edges. That is, the 
temperature rise profile over the panel surface will look like a truncated dome shape. High-intensity ther¬ 
mal loading could induce (1) thermal buckling, (2) material degradation, (3) thermal creep, (4) thermal 
yielding, (5) thermal cracking after cooling down, etc. Excess thermal deformation caused by thermal 
buckling could disturb the airflow field, creating localized hot spots that could degrade the panel's struc¬ 
tural performance. 
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When the structural panels are applied as the hypersonic aircraft wing skins, the aerodynamic loading 
during hypersonic flight will cause the wing upper panels to be under combined spanwise compression (re¬ 
sulting from wing bending) and shear (resulting from wing torsion). On the other hand, the wing lower skin 
panels will be subjected to combined spanwise tension and shear loading. Under thermal loading, both the 
wing upper and lower skin panels will be under mainly biaxial compression with certain localized shear. 
The thermal loading will increase the mechanical compressive stresses in the wing upper panels, and tend 
to reduce the mechanical compressive stresses in the wing lower panels. Thus, the thermomechanical 
buckling characteristics of the hot structural panels are a critical concern in the hypersonic aircraft wing 
structural panels. 

Thermal buckling problems of single plates (continuous or laminated composites) were investigated 
by several authors in recent years (refs. 1-6), and thermomechanical buckling characteristics of the hot 
structural sandwich panels were analyzed extensively by Ko and Jackson (refs. 7-12). Using the minimum 
potential energy method, Ko and Jackson developed thermomechanical buckling equations for orthotropic 
rectangular sandwich panels subjected to combined mechanical compressive and shear loading, or under 
thermal loading (refs. 11-12). 

This report investigates the thermal buckling characteristics of uniformly heated rectangular titanium 
sandwich panels of different aspect ratios supported under four different edge conditions. The thermal 
buckling loads (or buckling temperatures) will be calculated using Ko-Jackson thermal buckling equations 
(refs. 11-12) and the finite element method for the purpose of validating the Ko-Jackson theory. The ther¬ 
mal buckling solutions obtained from the two methods will be compared and their discrepancies discussed. 

DESCRIPTION OF PROBLEM 

Figure 1 shows the geometry of the hypersonic aircraft rectangular sandwich panel. The sandwich pan¬ 
el has length a, width b, and depth h . It is fabricated with titanium face sheets of the same thicknesses 
t s , joined together through a titanium honeycomb core of depth h c using the enhanced diffusion bonding 
process (ref. 13). Figure 2 shows the combined thermal forces acting on the middle plane of the sandwich 
panel. 

For the thermal buckling analysis, the panel will be subjected to uniform temperature field under four 
different edge conditions shown in figure 3. The edges in the x- and y-directions are defined, respectively, 
as sides and ends. The minimum potential energy thermal buckling theory developed by Ko and Jackson 
(refs. 11-12) and the finite element method will be used to calculate panel buckling temperatures, and the 
eigensolutions based on the two methods will be compared. 

RAYLEIGH-RITZ THERMAL BUCKLING ANALYSIS 

In the thermal buckling analysis of sandwich panels conducted by Ko and Jackson (refs. 11-12), the 
extensional and bending stiffnesses are provided by the two face sheets, and the panel transverse shear 
stiffness is provided by the sandwich core only. 
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Panel Boundary Conditions 

The four sets of boundary conditions used in the Ko and Jackson thermal buckling analysis 
(refs. 11-12) are given below. 

Case 1: Four edges simply supported (4S condition) 

x = 0, a: u = v = w = M x = y yz = 0 (1) 


y = 0, b: u = v = w = M y = y xz 

= 0 

(2) 

Case 2: Four edges clamped (4C condition) 



n <3w 

x = 0, a: u = v = w = — = v 
dx xz 

■ V ■ 0 

(3) 

n , dw 

y = 0, b: u = v - w - — = y 
dy * xz 

= V = 0 

(4) 

Case 3: Two sides clamped, two ends simply supported (2C2S condition) 


x = 0, a: u - v - w = M x = y yz 

= 0 

(5) 

n , dw 

y = 0, b: u = v = w = — = y 
dy 

II 

o 

(6) 

Case 4: Two sides simply supported, two ends clamped (2S2C condition) 


n dw 

x = {), a: u = v = w = — = y vr 
dx xz 

■ V - 0 

(7) 

y = 0, b: u - v - w - M y = y xz 

= 0 

(8) 


For anisotropic sandwich panels, cases 3 and 4 will give different thermal buckling solutions. 

Deformation Functions 

For satisfying the different sets of boundary conditions (1) through (8), the associated deformation 
functions { w, y xz , y yz } chosen by Ko and Jackson (refs. 11-12) in the thermal buckling analysis of the 
sandwich panels have the following forms: 

Case 1: 4S condition 


w (x, y) 




„ S 

b 


(9) 
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w(x,y) = sin^sin^ ^ J ( 12 ) 


(13) 


(14) 


(15) 


(16) 


(17) 
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Case 4: 2S2C condition 


, * . ;tx v, ^ . . rrmx . rmy 

w(x,y) = sin— > > A sin-sm—^ 

a x-i X-i mn a b 


(18) 


/ s TCX o mKX 

y^{ x >y) = cos — > > 5 sm—s 
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ju v r> nmx . rmy 
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, N . Jtx v, v, „ . rrmx 
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( 20 ) 


The choice of these four sets of deformation functions, each of which satisfies the associated boundary 
conditions (1) through (8), is for the mathematical amenability of the eigenvalue solutions. As shown in 
figure 4, the zero transverse shear distortions (i.e., y xz = 0 or y z = 0 ) at the panel edges cannot be en¬ 
forced simultaneously in the actual panel deformations, except for the 4C condition. 


Thermal Buckling Equations 


The thermal buckling equations developed by Ko and Jackson (refs. 11-12) for uniformly heated and 
constrained rectangular orthotropic sandwich panels using the Rayleigh-Ritz method, written in terms of 
temperature rise AT for each set of integral values {m, n} (or mode shape), have the following form 

| ^ [^f + P ’""» + = 0 (21) 

The bending-stiffness parameter M mnkl and the extensional stiffness parameter P mnkl in equation 
(21) are defined as 




-li 

a mnkl + 

classical 
thin plate 
theory term 


mnkS a mnkl a mnkl a rnnkfirnnkl 9 + a mnkl Ka mnkl a mnkl a mnkl Cl mnkl > 
22 33 23 32 

a mnkl a mnkl a mnkl a mnkl 


( 21 
Vfl 


J 


( 22 ) 


transverse shear effect terms 


P mnkl S ^ ab [S *) ^ A ll a x + A 12«v) + ^ Z) ( A ^ a x + A 22« ) ] (23) 

^66% 

The characteristic coefficients a ] mnkl (i, j = 1, 2,3) appearing in equation (22) are defined in appendix 
A, and a mnkl in equation (22) is the first part of a mnkl , containing no thermal loading terms (i.e., terms 
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containing {k x , k })• The parameters {£, £,} in equation (23) are, respectively, the numerical coeffi¬ 
cients of the load factors {k x , k y } contained in a^ M . The values of {£, £} change with the indicial and 
edge conditions (ref. 12). 

The numerical parameter r|, appearing in equations (22) and (23), and the special delta function 5 mM , 
appearing in equation (21), are defined for different edge conditions as (ref. 12) 

Case 1: 4S condition 


r\ = 32 

^mnki = —o— rnphL_— m ±k = odd, n±l = odd 


(24) 


Case 2: 4C condition 
(16) 3 

*1 = 


J mnkl 7 2 ' 9 9" ~ 2 9 2 2 ’ 

(m -k ) (n -l )[{m + k) - 4] [ (m - k) - 4] [ (« + /) - 4] [ (n -1) - 4] 


(25) 


m±k - odd, n±l = odd 


Case 3: 2C2S condition 
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^mnkl 


mnkl[2— {n + / 2 ) ] 


(m -k ) (n z -t)[(n + l)*-4][(n-iy-4] 


m±k = odd, n±l = odd 


(26) 


Case 4: 2S2C condition 
"r! = 8 3 




mnkl[ 2- (m + k ) ] 


(m -k)(n -l)[(m + k) -4][(m-k) -4] 


m±k = odd, n±l - odd 


(27) 


In the thermal buckling equation (21), both M mnk i and P mnk i terms contain temperature dependent 
material properties. Thus, in the eigenvalue solution process using equation (21), one has to assume 
a buckling temperature A T a and use the material properties corresponding to A T a as inputs to calculate 
the buckling temperature A T cr . This material property iteration process must continue until A T a 
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approaches A T cr Thus, in the thermal buckling, the eigenvalue solution process is a multi-step process. 
However, in the mechanical buckling, only one step eigenvalue solution process is required. 

When the coefficient of thermal shear distortion is zero (i.e., a xy = 0 ), the buckling equation (21) 
takes on the form 


?? 


M mnkl 

A T 


= 0 


where M mnkl and P mnkl 


M mnkl = 


a mnkl + 


are defined as 

12 , 23 31 

a mnkl ' a mnkl a mnkl 


21 33 . 13 ( 21 32 22 31 \“ 

~ a mnkflmnkl ' * a mnkl a mnkl a mnkl a mnkl a mnkl 
22 33 23 32 

a mnkl a mnkl a mnkl a mnkl 


(28) 


(29) 


Pmnkl *[l(m,k) (A n a x + A 12 a y ) + £ (n, l) (A 21 a x + A 22 a y ) ] (30) 

When equation (28) is reduced to the isotropic case with no transverse shear effects, the buckling tem¬ 
perature A T cr will be independent of the material’s modulus of elasticity (ref. 14). 

The characteristic equation (21) forms a system of infinite number of simultaneous homogeneous 
equations, each of which is associated with each indicial combination of {m, n}. Those simultaneous equa¬ 
tions, generated from equation (21), have the following characteristics. The first two terms 
(Mmnki !+ P m nkl ) °f equation (21) are nonzero only for the indicial conditions {m = k or m - k = 2} 
and {n = l or n -1 = 2} based on the indicial constraints for a J mnkl (appendix A). Thus, if (m ± n) is even, 
then (k ± t) is also even, and if (m ± n) is odd, then (k±l) is also odd. The special delta function § mnkl in 
the third term of equation (21) is nonzero only when (m ± k) is odd and (n ± 1) is odd. It follows that ( m ± k) 
± (n ± l) = (m ± n) ± (k ± l) = even. This implies that if (m ± n) is even, then (k ± t) is also even, and if 
(m ± n) is odd, so also is (k ± l). Because of these indicial characteristics, there is no coupling between the 
even case (i.e., symmetric buckling) and the odd case (i.e., antisymmetric buckling). Thus, the simulta¬ 
neous equations generated from equation (21) may be divided into two groups that are independent of each 
other—one group for which (m ± n) is even, and the other group for which (m ± n) is odd (refs. 11-12). 

For the deflection coefficients A kl to have nontrivial solutions for given aspect ratio b/a, the determi¬ 
nant of coefficients of unknown A^of the simultaneous homogeneous equations written out from 
equation (21) must vanish. The largest eigenvalue 1/A T thus obtained will give the lowest buckling tem¬ 
perature AT cr . The determinants of the coefficients of the simultaneous equations written out from equa¬ 
tion (21) up to order 12 are given in appendix B for the cases m±n = even (symmetric buckling) and m ± 
n = odd (antisymmetric buckling) for different edge conditions. The determinants of order 12 were found 
to give sufficiently accurate eigensolutions and, therefore, the determinants were truncated at order 12 in 
the present eigenvalue extractions. In appendix B one notices that for the 4S edge condition only, 
(Mmnkl! A^ + Pmnkl! terms form the diagonal terms of the determinants, and the nonzero off-diagonal 
terms consist only of numerical values given by & mnk! ■ However, for the rest of edge conditions, 
(M mnkl /A.T + P /nnk i) not only appear in the diagonal terms, but also in the off-diagonal terms (mixed with 
the numerical terms associated with & mnk[ )- 



FINITE ELEMENT THERMAL BUCKLING ANALYSIS 


The structural performance and resizing (SPAR) finite element computer program (ref. 15) was used 
in the finite element thermal buckling analysis of the sandwich panels. 

Finite Element Modeling 

To gather thermal buckling data of sandwich panels having a wide range of aspect ratios b/a, three ba¬ 
sic finite element models of different b/a were set up so that each model would cover certain limited range 
of b/a. Changing b/a of each model was accomplished by simply modifying length b and keeping length 
a constant. Sometimes more elements had to be added to an overstretched model to maintain proper ele¬ 
ment aspect ratios. An overstretched model without additional elements could result in local buckling of 
slender element cells rather than global buckling of the panels (i.e., local buckling temperature is less than 
the global buckling temperature). For low b/a (< 1.8) and high b/a (> 2.9) aspect ratio panels (figs. 5(a) 
and 5(c)) for which the lowest buckling mode were symmetric, only the quarter panels were modeled. The 
SPAR constraint commands SYMMETRY PLANE = 1 and SYMMETRY PLANE = 2 (ref. 15) were then 
used to generate the full panels. For moderate b/a (1.8 b/a s 2.9) aspect ratio panels (fig. 5(b)) for 
which the lowest buckling mode could be either symmetric or antisymmetric, half panels were modeled. 
The SPAR constraints command SYMMETRY PLANE = 1 was then used to generate the full panels. For 
purely antisymmetric buckling, one can model only a quarter panel and use the constraint commands 
SYMMETRY PLANE = 1 and ANTISYMMETRY PLANE = 2 to generate the whole panel. However, 
such quarter-panel models were not used in gathering the buckling data because they consistently gave 
somewhat higher buckling temperatures than those given by the half-panel models. Figure 6 shows the 
three basic finite element models set up for the sandwich panels. Both models A and C are the quarter- 
panel models, but model B is a half-panel model. From these basic models, several modified models (not 
model shown) were also set up for handling certain aspect ratios and edge conditions. 

Each face sheet of the sandwich panel was modeled with one layer of E43 elements (quadrilateral com¬ 
bined membrane and bending elements) and the sandwich core with one layer of S81 elements (hexahe¬ 
dron (or brick) elements), which connect the upper and the lower face-sheet elements E43. Because the 
joint locations of those face sheet elements E43 are located in the middle planes of the respective face 
sheets, the finite element core depth will then be h instead of the actual depth h c (fig. 7). Thus, to simulate 
the actual relative displacements (or maintain same stiffness) between the two face sheets in the sandwich 
thickness direction and the x- and y-directions (i.e., d 3 = d' 3 and d x = d\ and d 2 = d' 2 , fig. 7), the 
thickness elastic modulus E' cz , and the transverse shear moduliG^. 2 and G' cyz of the S81 elements had to 
be increased slightly according to the following relationships: 


E'cz 

II 

(31) 

C'cxz 

II 

(32) 

G cyz 

r h 
= a ‘y>h 

(33) 
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One can also model the sandwich core with one layer of S81 elements having the exact depth h c , and 
then connect the gaps between E43 joints and S81 joints with rigid elements. However, this alternative 
modeling method requires twice as many total joint locations, and therefore, it was not used. 

For simply supported edges, free rotation and free transverse shear deformation must be allowed 
(fig. 8(a)). To simulate this type of edge, pin-ended rigid rods were attached to the panel edge for connect¬ 
ing the two face sheets (sandwich core carries no extensional and bending stiffnesses), and then the mid¬ 
points of the rigid rods were pin-jointed to fix points lying in the sandwich middle plane (fig. 8(a)). Each 
pin-ended rigid rod was modeled with two identical E22 elements (beam element for which the intrinsic 
stiffness matrix is given). To simulate the rigidity of the rod, extensional and transverse shear stiffnesses 
of the E22 elements were made very large. The pin-joint condition at the face sheet edges was simulated 
by assigning zero values to the rotational spring constants in the stiffness matrix for the E22 elements. The 
pin-joint condition at the middle-plane fixed points was simulated by relaxing the three rotational con¬ 
straints. Two methods were used to connect the ends of the E22 elements to the middle plane fixed points. 
In the first method (center drawing of fig. 8(a)), the first joint of each E22 element was connected to the 
associated joint of E43 element and its end point to the panel middle-plane fixed point. In the second meth¬ 
od (right-hand drawing of fig. 8(a)), the ends of the two identical E22 elements, whose first joints were 
connected to the upper and the lower face sheets, were connected together to the middle-plane fixed point 
through E25 element (zero length element used to elastically connect geometrically coincidental joints). 
The stiffnesses of the E25 were made so large that the two E22 elements, connected together by the E25 
element, will behave like one rigid rod. These two types of simply supported edge simulations were found 
to give identical thermal buckling solutions. In most of the buckling data gathering, the first edge simula¬ 
tion was used because it required less joint locations. 

For the clamped edge (fig. 8(b)), the edges of the two face sheets were built into fixed vertical walls to 
generate the desired constraints of zero slope, zero in-plane displacements, and zero transverse shear de¬ 
formations. 

Table 1 shows the sizes of the three finite element models set up for the sandwich panels of different 

b/a. 


Table 1. Sizes of three finite element models A, B, and 
C (c.f., fig. 6). _ 


Feature 

Model A 

Model B 

Model C 

JLOC 

2178 

4850 

4050 

E43 

2048 

4608 

3840 

S81 

1024 

2304 

1920 


To make sure that the above finite element models gave accurate eigensolutions, the sandwich cores 
of model A was modified to two layers of S81 elements to investigate the eigensolution convergencies. It 
turned out that both the basic and the modified models gave practically identical eigensolutions. Because 
the modified model required about three times longer computational time, it was not used in the actual 
buckling data gatherings. 


11 



Eigenvalue Extractions 

The eigenvalue equation for buckling problems is of the form 


\KX + KX = 0 (34) 

where 

K g = system initial stress stiffness matrix (or differential stiffness matrix), corresponding to partic¬ 
ular applied force condition (e.g., thermal loading), and in general a function of X 

K = system stiffness matrix 
X = displacement vector 

= eigenvalues for various buckling modes 

The eigenvalues X.. (i = 1,2,3,...) are the load factors by which the static load (mechanical or thermal) 
must be multiplied to produce buckling loads corresponding to various buckling modes. Namely, if the ap¬ 
plied temperature load is AT, then the buckling temperature A T cr for the i'-th buckling mode is obtained 
from 


A T cr = X.A T (35) 

If it is desired to find eigenvalues in the neighborhood of c, then the following shifted eigenvalue equation 
may be used. 


( X-c)K g X+ (K + c K g )X = 0 (36) 

In the eigenvalue extractions, the SPAR program uses an iterative process consisting of a Stodola ma¬ 
trix iteration procedure, followed by a Rayleigh-Ritz procedure, and then followed by a second Stodola 
procedure. This process results in successively refined approximations of i eigenvectors associated with 
the i eigenvalues of equation (34) closest to zero. Reference 15 describes the details of this process. 

NUMERICAL EXAMPLES 

The hot structural sandwich panels analyzed were fabricated with titanium face sheets and titanium 
sandwich core having the following geometrical and material properties. 

Geometry: 


a =24 in. 
b = varying 
h = 0.75 in. 
h c = h-t s = 0.69 in. 
t s = 0.06 in. 
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Material properties: 


Face sheets 



70 °F 

900 °F* 

E x = E y , lb/in 2 

16 x 10 6 

13.1 x 10 6 

lb/in 2 

6.2 x 10 6 

5.0 x 10 6 

V = V 

xy yx 

0.31 

0.31 

a x = a y , in/in-°F 

4.85 x 10 -6 

5.35 x 10~ 6 

a xy , in/in-°F 

0 

0 

p Ti , lb/in 3 

0.16 

0.16 


Mach 15 flight temperature. 


Honeycomb core (properties at 600°F) 



2.7778 x 10 4 lb/in 2 

2.7778 x 10 4 lb/in 2 

2.7778 x 10 5 lb/in 2 
0.00613 lb/in 2 
0.81967 x 10 5 lb/in 2 
1.81 x 10 5 lb/in 2 
0.658 x 10~ 2 
0.643 x 10 -6 
0.643 x 10 -6 
5.37 x 10" 6 in/in-°F 
0 in/in-°F 
3.674 x 10- 3 lb/in 3 


The main objective of the present report is to study the general trend of thermal buckling characteristics 
of sandwich panels under different edge conditions, and to validate Ko-Jackson theory (ref. 12). Because 
of the lack of material property data at high temperatures, the material iteration process was not performed, 
and the face sheet properties at 900 °F and sandwich core properties at 600 °F were used in the buckling 
temperature calculations. 


13 



RESULTS 


Eigenvalue Iterations 

In finite element eigenvalue extractions, the maximum number of iterations was set to be 100. For most 
cases (with or without eigenshifting), however, the eigenvalues converged well below 100 iterations based 
on the convergence criterion (| X,;- X,/_i | / X,,) < 10~ 4 . Figure 9 compares the convergence curves of eigen¬ 
value iterations with and without shifting for the square panel (b/a = 1) under the 4C condition. With shift¬ 
ing, the number of iterations could be reduced from 14 to 9 iterations. For certain cases, the number of 
eigenvalue iterations with shifting turned out to be very close or even identical to that without shifting. For 
certain problems, such as the thermocryogenic buckling of cryogenic tanks (ref. 16), eigenshifting could 
drastically reduce the number of eigenvalue iterations (i.e., reduction in computer time). For the present 
sandwich panel buckling problem, however, the reduction in the number of eigenvalue iterations through 
the eigenshifting turned out to be relatively small or negligible. 

In most of the thermal buckling data gathering, the eigenshifting method was used. The shifting factors 
used were near the values of the buckling temperatures predicted from the energy theory. In figure 10 the 
ELXSI 6400 computer processor times are plotted as functions of the number of eigenvalue iterations for 
the four edge conditions. The processor time per iteration (i.e., slope of the data fitting line) for the 4S, 4C, 
2C2S, and 2S2C conditions are, respectively, 2.45,1.88, 2.48, and 2.06 minutes. 

Buckling Temperatures 

Figures 11,12, and 13, respectively, show the lowest-mode buckling shapes of sandwich panels of as¬ 
pect ratios b/a = 1, 2, and 3 under the four different edge conditions. Notice that at the simply supported 
edges of the 4S case and at the clamped edges of the 2C2S and 2S2C cases, the transverse shear deforma¬ 
tions cannot be zero (fig. 4). At higher b/a (figs. 12 and 13), the 4C and 2S2C cases required finer element 
models for obtaining global buckling. The square panel (b/a = 1, fig. 11), under all the four different edge 
conditions, buckled symmetrically with {m = 1, n = 1} buckling mode. For the b/a = 2 rectangular panel 
(fig. 12), the 4S, 4C and 2C2S conditions still induced symmetrical buckling mode of {m = 1, n = 1}. How¬ 
ever, under the 2S2C condition, the rectangular panel buckled antisymmetrically under {m = 1, n = 2} 
buckling mode. For the slender panel of b/a = 3 (fig. 13), both 4S and 2C2S conditions continued to induce 
symmetrical buckling mode of {m = 1. n = 1}. However, under the 4C and 2S2C conditions, the multiple 
symmetrical buckling mode of {m = 1, n = 3} turned out to be the lowest buckling mode. In figure 14 the 
thermal buckling temperatures calculated using the minimum energy method (solid curves) and the finite 
element method (circular symbols) are plotted as functions of panel aspect ratio b/a for the four cases of 
edge conditions. Notice that for the high values of b/a, the thermal buckling solutions for the 2C2S and 
2S2C cases approach those of the 4S and 4C cases, respectively, because the constraint effects of the short¬ 
er edges of the slender panels diminish. The buckling solutions obtained from the two methods compare 
fairly well. The average solution difference between the two methods are 3.87%, 1.62%, 2.04% and 2.71% 
respectively for the 4S, 4C, 2C2S and 2S2C cases. The finite element method tends to give slightly lower 
buckling temperatures than those given by the minimum-energy method. The reason could be the follow¬ 
ing: (1) the finite element method allows deformations in the sandwich thickness direction, which the min¬ 
imum energy theory ignores, (2) the theoretical edge conditions of zero y xz and y yz (eqs. (1), (2), (5)-(8)) 
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cannot be enforced properly in the finite element edge constraints for the 4S, 2C2S, and 2S2C conditions 
(fig 4), and (3) the finite element modeling assumptions. For the 4C case only, all the theoretical edge con¬ 
ditions could be enforced in the finite element edge constraints. For 4S, 4C and 2C2S cases, the discrep¬ 
ancy of the eignsolutions between the two methods is larger at the low values of b/a, and gradually 
diminishes at high values of b/a. This solution discrepancy is minimum for the 4C case, and maximum for 
the 4S case (because both y xz and y z cannot be zero at the edges of the finite element modes). For the 
2S2C case, the solution discrepancy is almost unaffected by the change of b/a. Table 2 lists the buckling 
temperatures A T cr calculated from the minimum energy and the finite element methods. 


Table 2. Buckling temperatures of sandwich panels calculated using minimum 
energy and finite element models. _ 


— 




A T cr , 

°F 






4S 


4C 

2C2S 

2S2C 

b/a 

Energy 

theory 

Finite 

element 

Energy 

theory 

Finite 

element 

Energy 

theory 

Finite 

element 

Energy 

theory 

Finite 

element 

0.5 

1297 

1207 

2541 

2498 

2424 

2366 

1512 

1464 

0.6 

1051 

970 

2169 

2126 

1995 

1938 

1334 

1286 

0.7 

885 

815 

1889 

1847 

1654 

1603 

1232 

1187 

0.8 

769 

710 

1683 

1645 

1387 

1343 

1175 

1132 

0.9 

682 

637 

1535 

1499 

1181 

1143 

1144 

1105 

1.0 

622 

583 

1428 

1396 

1021 

988 

1128 

1093 

1.1 

575 

547 

1352 

1322 

897 

866 

1122 

1090 

1.2 

538 

513 

1297 

1271 

799 

774 

1121 

1092 

1.4 

486 

465 

1232 

1209 

662 

640 

1126 

1102 

1.6 

451 

436 

1199 

1179 

573 

560 

1136 

1116 

1.8 

427 

416 

1183 

1165 

514 

504 

1144 

1130 

2.0 

409 

403 

1175 

1160 

473 

471 

1128 

1096 

2.2 

396 

391 

1172 

1158 

444 

439 

1122 

1092 

2.4 

386 

382 

1171 

1156 

423 

419 

1121 

1094 

2.6 

379 

376 

1164 

1155 

407 

404 

1122 

1099 

2.8 

372 

370 

1156 

1151 

395 

393 

1126 

1100 

3.0 

367 

366 

1150 

1145 

386 

384 

1128 

1095 

4.0 

353 

348 

1140 

1121 

360 

360 

1124 

1100 


CONCLUSIONS 

Thermal buckling characteristics of hypersonic aircraft honeycomb-core sandwich panels subjected to 
uniform temperature loading were analyzed using minimum energy theory and finite element methods. 
The thermal buckling curves were generated for titanium sandwich panels of various aspect ratios. The 
two methods predicted very close buckling temperatures, and thus, the Ko-Jackson theory was validated. 
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The finite element method tended to give slightly lower buckling temperatures than those given by the 
minimum energy theory. The slight discrepancies in the eigensolutions between the two methods could be 
attributed to the following: 

1. The minimum energy theory does not consider deformations in the panel thickness direction, where¬ 
as the finite element method does. 

2. The theoretical zero transverse shear deformations at the panel edges cannot be enforced in the finite 
element models with simply supported edges and cannot be enforced simultaneously in the finite elements 
models with mixed simply supported and clamped edges. 

3. Assumptions made in finite element modeling. 

The discrepancy of the eigensolutions between the minimum energy theory and the finite element 
method is the largest for the simply supported edge condition, because the zero transverse shear deforma¬ 
tions at the panel edges cannot be constrained in the finite element models. This solution discrepancy is 
larger at low values of b/a and gradually decreases as b/a increases. For the sandwich panels the eigen- 
shifting has small effect on the improvement of the eigenvalue convergence rate. 

The author gratefully acknowledges the contribution of the late Raymond H. Jackson, NASA mathe¬ 
matician, in setting up computer programs for the eigenvalue extractions. 

Dryden Flight Research Center 

National Aeronautics and Space Administration 

Edwards, California, June 14, 1994 
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940463 


y 


(a) Four edges simply supported (4S condition). 


y 940464 

(b) Four edges clamped (4C condition). 



(c) Two sides clamped, two ends simply support¬ 
ed (2C2S condition). 



Figure 3. Four types of edge conditions. 
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(c) Model C—Quarter-panel model for {m = 1, n = 3}. 

Figure 6. Three finite element models generated for sandwich panels of different aspect ratios 





Actual panel 


3 = d 3 —|- 



d i = d i 

d 2 = d 2 — 


Figure 7. 


Finite element model 



of sandwich panel. 



(a) Simply supported edge. 



Figure 8. Simulation of different edge conditions. 
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Figure 9. Convergence curves of eigenvalue iterations; 4C condition; b/a = 1. 
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Figure 10. Increase of processor time with number of eigenvalue iterations; ELXSI6400 computer. 
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Figure 14. Buckling temperature curves for titanium sandwich panels under different edge conditions; 
a = constant. 




APPENDIX A 


COEFFICIENTS OF CHARACTERISTIC EQUATIONS 


The characteristic coefficients a l, mnkl appearing in equation (22) are defined in the following for differ¬ 
ent indicial and edge conditions (ref. 12). 

Case 1: 4S condition 


(1) m = k, n = l 



(2) m k k ,n k l 


(A-l) 


a mnkl 


= 0 


(A-2) 
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Case 2: 4C condition 


11 

a llll : 

= 12D .i(f) 

4 + 8 (£>12 

^«)(S 2 (S 2 

+UI> »(s ) 4 



~M=) 

l W=)’ 
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a llll : 

21 

= "mi ~ “ 

[ 12D .*(; 

;) 3 + 4(C, 2 + 2£) 66 : 

>(;)(s) 2 ] 
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a llll : 

31 

= "mi = “ 

[l 2 Ml 

;) 3 + 4(fl 12 + 2i) (i6 ; 
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(A-3) 

22 

"llll : 

■ 12D *0 

2 + «>«( 

I) 2+3 ^ 




■m, =“m. =“(D 12 + fl e6 )(5)(5) 
•mi = UDj^\iDj$\lD Qy 


(2) m = k = 1, n = l V 1 


, = 8£>„(?) + 4(O i2 + 2D 66 )(5) (5) (1 + „ 2 ) + |b 22 (?) [ (1 + n 2 ) 2 + 4n 2 ] 

, = 41l. = - [ 8i> .l( f ) 3 + 2 + 2D ^) ( f) ( S) 2< 1 + " 2) ] 

. = “ill, = - ||fl 2 2 ( l) 3 [ ('1 + "V 2 + 4n 2 ] + 2 (D 12 + 2 DJ ( f) "( j) (1 + n) J 
, = 8D, 1 (5) 2 + 2C 66 (!) 2 (l +n 2 ) + 2B & 


33 3 „ / Jt' 
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(3) m = k k 1, n = l = 1 


|c„(5) 4 [(1 + m 2 ) 2 + 4m 2 ] + 4 (D 12 + 2Z> 66 ) (5) 2 (5) 2 < 1 + m 2 ) t 8^45) 4 

= -||o„(j) 3 [(l + m 2 ) 2 + 4m 2 ] +2(0 12 + 2Z> S6 )(5)(5) 2 (l+m 2 )J 

“l\.l = - [8022( I) 3 + 2 (f 12 + 2D m ) ( 5) 12 ( j?) (1 4 m 2 )] 

l-Oi. (f) :2 [ ( 1 + •"S ' 2 + 4 » 2 ] + 2«4 1) 2 <1 4 ™ 2 ) + |o e , < 1 + m 1 ) 

4L, = 2<2>,2 + »66)(;)(5)(l+'» 2 ) 
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(5) m = k = 1, n -1 = 2 


“Wu = - [«>„(;)'* + 2<o 12 + 20«) (f )*( s) 2 <' 1 + 0 :2 + - 4 d 22 (!)' 4 (:l + /)'‘] 

+fl '(=) 2 [*a 2+ Hf) 2 < 1+ '> 2 ] 

= «!',./ = iMf) 3(1 + ,)4+ (^ 12 + 2r. 6)i )(=) 2 (|) (1+/) 2 (A-7) 

“llu = - [ 4£> n(S) Z+jt> «(i) 2 (i + 0 2 +i> e J 

= -< c i2 + °66>(;;)(£)<i + '> 2 

- [iM s) 2 < 1 + ()4+c 4 j) 2 < 1 + ()2 + 5 fl e, < 1 + ') 2 ] 

(6) m = k k 1, n-l = 2 

«1w = 4{°4S 4 [< 1+ '" 2 > 2 + 4 ” 2 1 +2 ®u* 2 o«)(2) 2 (I) 2 (1 + ». 2 )(i + 0 2 
+ D45) 4 ( 1 + /) 4 }4(f) 2 [^(5) 2 (l + m 2) + ^(I) 2 ( 1 + 0 2 ] 

“W, - ‘1., - ||o 11 (j) 3 [<1 + "' 2 ) 2 + 4». 2 ] ♦ (O 12 + 20 66 )(5)(l) 2 d +m 2 ) <l + /) 2 j 

= “»,»< = !{°4l) 3(i+ '> 4 +(d >2 +2b «.> 0 2 (5) <1+m2) (1+ ' )2 } (a-8) 

«™-, = - 5 { u ( J) 2 [ l1 + n ) 2 + ] + M J) :2 c1 + '» 2 ) (■1 + '>' 2 + D Qx (1 + m 2 )} 

«»„/ = «L,/ = (0 12 +B 66 )(5)(l)(l + m 2 )(l + ;) 2 

= -i [°» (s) 2(1 + i)4+D « (^) 2(1 + m2) ( 1 + ; ) 2 +o s ,(i + 0 : 


35 



(7) m-k = 2, n = l = 1 
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(9) m-k = 2, n-l = 2 
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(2) m = k, n = l k 1 
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Case 4: 2S2C condition 
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(3) m-k = 2, n = l 
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APPENDIX B 


BUCKLING EQUATIONS 

The buckling equations (eigenvalue solution equations) written out from equation (21) up to order 12 
(i.e., 12 x 12 matrices) for the cases m±n = even (symmetric buckling) and m±n = odd (antisymmetric 
buckling) for different edge conditions are given on the following pages (ref. 12). 
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= odd (antisymmetric buckling) (4S) 







Case 2: 4C condition 
m ± n = even (symmetric buckling) 




~W 1315 4725 


~~AT~ r 1333 33075 


_ 2024 M 335 

33075 Ar 


AT * 2444 H025 


i=5, n =3 



= odd (antisymmetric buckling) (4C) 




Case 3: 2C2S condition 
m ± n = even (symmetric buckling) 




= odd (antisymmetric buckling) (2C2S) 





Case 4: 2S2C condition 
m ± n = even (symmetric buckling) 




= odd (antisymmetric buckling) (2S2C) 
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